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ABSTRACT 


For  any  n-times  differentiable  function  f with  uniform  bounds  on  f and 
f^n\  we  study  the  pair  of  values  (f^(t),  f^  ^ (t)  ) for  an  arbitrary  real 
t and  a prescribed  j = 0,...,n  - 1.  A given  value  of  f^(t)  determines 

(j  + 1) 

admissible  values  for  f (t) . These  values  are  exactly  determined  in  terms 

of  the  Euler  spline  (t).  Special  differentiation  formulas  of  cardinal  inter- 

n 

polation  type  are  developed  to  solve  the  problem. 
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SIGNIFICANCE  AND  EXPLANATION 


A simple  physical  application  of  the  problem  solved  in  this  paper  is  the 
following:  Given  a motion  f(t)  constrained  to  lie  always  in  a certain 
interval  and  always  with  a bounded  acceleration  f"(t),  best  possible  bounds 
for  the  velocity  f' (t)  can  be  given.  It  is  intuitively  clear  that  the  bound 
on  the  velocity  f'(t)  at  time  t depends  on  the  position  f(t).  We  study 
this  situation,  and  more  generally  the  case  of  pairs  of  derivatives  of  higher 
order. 
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A RFldNLMFNT  t i KOI MilGOROV'S  inequality 


A.  S.  Cavaretta,  Jr. 


§ 1 . Inti  t »Juction . 

In  1939  Kolmogorov  [4]  proved  a sharp  inequality  between  the  suj  remum  norror  of  the 
successive  derivative  of  a function.  With  n >_  2 and  valuer-  for  |j  f |J  and  ||f  (n)  || 
he  found  best  possible  estimates  for  |,1  ^ j|  , 1 5 j 1.  n ~ ^ * here,  and  in  all  that 
follows,  tin*  norm  r s the  supremum  norm  taken  over  the  entire  real  axis.  The  inequality 
is  intrinsically  ti'd  up  with  the  > called  Filler  ; pi  i ne  function  F^  (s)  and  can  be 
considered  as  a charact.  1 i 1 ic  j roper ty  of  F ^ - In  fact,  if  we  set 

- II  I! 


jn 


j = 1 n 


then  the  Kolmogorov  Theorem  tala,  or:  tire  following  form: 

Support  f h o in  a*-ol  nt  el  v ec-nt.j  nuous  (n-1) th  d_er  iv.pt  iyc  and  rati  r.f  ie: 


(1.1) 


I'l!'  !'f  I!  - Y,m  • 


Then  also 

U.2)  Ilf'1’  II  Yjn-  3 1 n - ) . 

Then*  i 1. i t i < , • 1 • ? ; ibS  as  they  are  < foi  F^(  s)  . 

The  constants  Y tan  1 r « idily  computed  from  the  Fourier  series  of  F (s)  . l\ 
jn  n 

a change  cf  s 1 1 * in  loth  tx*  • can  always  arrange  to  have  (1.1)  for  any  given 

function  f.  The  Filler  spline  occurs  most  naturally  within  the  context  of  cardii..»i 

n 

spline  i nt.orpolat  i n v*;.-.  :e  it  appears  as  the  unique  interpolant  of  the  sequence  ( - 1 ) • 

wc*  refci  tile  reader  to  I't)  page-:  39-40  and  also  to  [81  for  background  information  on 

these  remarkable  functions.  Most  importantly,  we  need  the  following  three  properties 

of  F (s)  : 
n 

i)  ||ffn(s)||  1 

(1.3)  ii)  fc’  (\>)  ( - 1 ) V for  all  inlpgois  v; 

n 

iii)  if  n is  oven  and  v any  integer. 
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when  v 


- + v 

(-1)  2 fc’‘n>  (s)  = ■> 

n nn 


1 1 
" 2 < S * V + 2 


if  n is  odd  and  v any  inteqer, 
n-1  . 

~ + v . . 

(-1)  E 11  (s)  = y when  v - 1 < s < v . 

n nn 


These  properties  characterize  i?  (s)  and  are  sufficient  for  our  purposes. 

For  convenience,  let  us  denote  by  F all  those  functions  f satisfying  the 

n 

hypotheses  (1.1)  of  the  Theorem.  Now  for  each  j = 0,1,... ,n  - 1,  define 

A.  = { (f  (j)  (s)  , f ( j + 1)  (s) ) ) 

3 

where  f ranges  over  the  whole  class  F and  s ranges  over  the  whole  real  axis. 

n 

Since  F is  invariant  under  shifts  of  origin,  we  may  set  s = 0 or  any  prescribed 
n 

value  t if  convenient.  When  we  view  A as  a subset  of  the  x - y plane  with 

,(j)  . . , , ,(j+l)  , , 

x = f (s)  and  y=f  (s)  , 

several  geometric  features  become  immediately  obvious.  Each  A is  convex.  Also  as 

f t F implies  if (is)  f F , we  easily  establish  that  A.  is  symmetric  in  each  axis, 
n n j 

And  from  the  Kolmogorov  Theorem  we  conclude  that  A is  a bounded  set;  more  precisely, 

it  is  circumscribed  by  the  rectangle  determined  by  the  lines  x = iv.  and  y = ±v . 

]n  * '3+1, n 

A comp  letc^de. script  ion  of  A is  given  by  the  following 


Theorem 


1 : I.et  0 < j < n - 2.  Th(‘  boundary  of  A is  given  pa r amc trie-ally  i n t by 


the  curve 


x ( t)  --  K<:i)  (t) 
n 

y ( t)  = F (:j41)  (t)  . 
n 


Since  E (t)  is  periodic  with  period  2 the  boundary  of  A is  parameterized  over  the 
finite  interval  [ 0 , 2 J and  is,  of  course,  a simple  closed  curve.  For  j = 0,  the 
result  is  already  implicit  in  Kolmogorov's  paper  of  1939  (4).  This  case  is  formulated 
there  as  an  auxiliary  inequality  used  in  the  induction  proof  of  the  tain  result  (1.2) 


-2- 


on  norm  inequalities. 


The  case  j - n - 1 is  exceptional  in  that  A 


reduce;;  to 


a rectangle. 

j = 1 

use  these  to 


n-1 


The  contribution  of  the  present  paper  lies  in  its  methods  and  the  cases 


2.  In  §2  wo  present  certain  interpolation  formulas  of  cardinal  type  and 
give  a proof  of  Theorem  1.  We  derive  these  formulas  in  §3. 


52. 


Some  formulas  of  cardinal  type;  a proof  of  Theorem  1 . 


We  could  define  the  sets  A.  for  function  classes  other  than  F . For  example, 

1 n 

let  B denote  all  entire  functions  of  exponential  type  n which  when  restricted  to 

IT 

the  real  axis  are  uniformly  bounded  by  1.  As  above,  put 

A = { (f (s) ,f * (s) ) |f  e , s real)  . 

For  A we  have  a 

Proposition.  The  boundary  of  / ££  given  parametrically  by  the  curve  (cos  n(  -n  sin  7T t ) . 

This  proposition  is  implicit  in  earlier  work  of  Duff  in  and  Schaeffer  (3),  and  indeed 

follows  quite  easily  from  a formula  of  Polya-Szego  (7;  III,  165].  Our  use  of  this 

formula  demonstrates  the  method  by  which  we  will  derive  our  Theorem  1. 

Proof.  We  exploit  the  following  formula,  valid  for  any  f e B and  any  t,  real  or 
n 

complex: 

OU  2 

(2.1)  TTCOSTTt  f (t)  - sinTTt  f * (t)  = y (-1)  --1--  f (v) 

» £ 

v=-OD  (t-v) 

CO 

= l A f(v) 
v=-«  v 

where  the  last  equality  merely  server,  to  define  the  coefficients  of  the  formula. 

Note  that  when  t is  real 

v 

siqn  A = (-1) 
v 

unless  t is  an  integer  for  then  all  but  one  of  the  A vanish. 

v 

Now  as  in  the  introduction  A is  viewed  as  a convex  subset  of  the  x - y plane. 

So  A has  a supporting  line  with  normal  vector  (a, 8),  see  Figure  1,  and  the  position 
of  this  line  is  determined  by 

(2.2)  max(ax  + 8y  [(x,y)  e A)  . 

Setting 

(2.3)  a = n cosii  t , 8 = -sinut 

for  an  appropriate  t,  we  see  that  the  corresponding  quantity  in  (2.2)  becomes 

(2.4)  rcosirt  f (s)  - sinnt  f'(s) 


-4- 


y = f 1 (s) 


- X = f (s) 

(a , R) 


Figure  1 

which  must  be  maximized  over  all  f e B and  over  all  real  s.  But  B is  invariant 

IT  71 

under  shifts,  so  we  may  just  as  well  take  s = t in  (2.4)  and  so  recover  the  left  hand 
side  of  (2.1).  Given  the  alternating  signs  of  A^,  formula  (2.1)  then  makes  clear  that 
(2.4)  with  s replaced  by  t is  maximized  when  the  function  f(s)  is  costts;  hence 

(cOSTTt,  -7TSin7Tt)  e 3/1  . 

This  result  persists  for  every  t,  and  varying  t we  generate  every  normal  direction 
(a,B)  as  seen  from  (2.3).  Thus  (cositt,  -ns inn t)  describes  the  full  boundary  of  A, 
as  was  to  be  shown. 

After  this  short  digression,  we  return  to  our  main  interest:  the  function  class  F 
and  the  corresponding  sets  A.,  j = 0,  , . . , n - 2 . Our  main  goal  is  a class  of  formulas 
analogous  to  (2.1).  The  existence  and  character  of  those  formulas  is  the  content  of 
Theorem  2. : Fix  n and  j with  n > 4 and  0 < j < n - 2 . Also  f i x a real  value  t. 


Then 

for  ai 

(2.5) 

n 

where 

(t)  = l A f (v)  + / K(s)f  (n)  (s)ds 


i)  (-1)  A > 0,- 

V 

ii)  K(s)  is , except  for  a discontinuity  at  t,  a cardinal  spline  of  degree 
n - 1 with  knots  a_t_  the  integers;  the  discontinuity  at  t is  in 
K(n'j'1)  and  K(n-j'2); 


iii)  for  n even 


-5- 


v+  — 

(-1)  2 K ( s ) > 0 i_f  v-^<s<v4- 


for  n odd 


V+ 


n-1 


(-1)  K (s)  > 0 i_f  v - 1 < s < v ; 

iv)  both  and  K(s)  tend  exponential ly  to  0 as 

tend  to  infinity. 


and 


Remarks.  The  and  K(s)  both  depend  of  course  on  t.#  but  we  do  not  indicate  this 


in  the  notation.  Formula  (2.5)  is  valid  for  every  f with  f 


(n) 


essentially  bounded 


and  f^n  ^ absolutely  continuous.  The  case  j = 0 and  t.  an  integer  is  exceptional 
as  then  the  left  hand  side  of  (2.5)  collapses  to  a multiple  of  f(t). 

The  proof  of  Theorem  2,  which  is  technically  compl icated , we  defer  to  53.  Here 
instead  we  give  in  detail  some  special  cases  and  then  indicate  how  (2.5)  and  (2.G)  are 
u‘  prove  Theorem  1.  We  observe  that  the  very  existence  of  formula  (2.5)  with 

•s  (2.6)  is  enough  to  establish  the  extremal  property  of  E^(s)  given  in 


For  our  first  example  of  the  type  of  formula  contained  in  Theorem  2,  set  n = 3 
and  j=0.  We  find  that  for  0 < t < 1 


(2.7) 


£2<t+l)f,(t)  " £’2(t+  2 >f  (t)  = 4(t-l)2€<0) 


1 


4t/f  (1)  4 / K (s)f  (3)  (s)ds 
0 1 


where 


Kt<s) 


? 2 

2<  t-irs 


2t2(S-n2 


0 < s < t 


t < s < 1 


Note  tha 1 K^ts)  >_  0.  There  arc  formulas  similar  to  (2.7)  for  other  values  of  t;  but 
due  to  the  symnetr i es  of  , (2.7)  is  sufficient  for  our  needs. 

For  n = 3,  j = 1 and  0 < t < 1 the  required  formula  is 


(2.R) 


£I(t+bf"(t)  - £H(t+i)f  ' «t)  = 8f  (0)  - Of  (1 ) 4 / K (s)  f (^)  (s)  ds 
2 2 7 2 „ t 


where 


-6- 


4s 


0 < s < t 


Kt(s)  = 


4 (s-1)  t < s < ] 

When  t = 0,  wo  infer  by  continuity  that  the  coefficient  of  f'(t)  in  (2.8)  is  -8; 
(2.8)  thus  icauces  to  the  Taylor  expansion  for  f ( 1 ) about  the  origin. 

The  case  n = 3 and  also  n = 2,  which  we  omit,  are  exceptional  in  that  our 
formulas  are  finite  in  nature.  The  situation  changes  for  n > 4,  for  then  we  have  the 
full  force  of  Theorem  2 and  the  formulas  arc  truly  of  cardinal  type,  involving  all 
integers  v as  nodes  and  kernels  K(s)  supported  on  the  entire  real  axis.  The  first 
sue)]  we  encounter  is  for  n = 4 and  j = 0: 


(2.9) 


OO  oo 

£3<t+2>f'(t)  ' £’i(t+2)f<t)  = I A f<v)  + / K(s)  f (4)  (s)ds 


i ) A. 


fll  (t)  ' v i 1 


a2ttU2' 


v < -1 


A,  = -11  + 2/30  ---  -.045548,  A _ = A 
1 2 1 


-1 


(2.10) 


ii)  for  0 < t • - and  M 


( 1 — A ) 


1 4 A, 


a . (t ) 


A 1 ' 1 2 2 

- t (4 1 13)  + ( 1 ) 

1 


. 1 4 A. 


1 8t3 


1,2 


A = A (t)  ■ -3  (4 1 2— 1 ) + yt2(4t2+3) 


(-1)  K(s)  > 0 for  V - 
An  easy'  calculation  from  ii)  shows  that 


1 1 

— < s < v + — . 
2 2 


a . (t)  > 0,  i = 1,2  and  A (t)  > 0 ; 
l 0 


hence  with  > 


0,  i = 1,2,  i)  implies 


(-1)  A >0  for  all  v 
v 
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1 

I 


Concerning  the  sign  regularity  of  K(s)  given  by  iii),  we  make  a series  of  remarks. 
From  our  construction  of  K(s)  in  §3,  it  will  be  clear  that  K(s)  has  simple  zeros 


at  every  point  v + — 
parts  the  remainder 


It  follows  that 


and  in  particular 


Once  K is  constructed  (2.9)  emerges  when  we  integrate  by 


/ K(s)f  <4)  (s)ds  . 


K'  ' ■ (1+)  - K'  ' ' (1-)  = -Aj  >0 


K"  ' (s)  > 0,  1 < s < 2 


K'"m  > 0 ■ 


Anticipating  considerations  of  §3,  we  find  three  (weak)  sign  changes  in  the  sequence 


K"  - , K' 


which  then  forces 


K'U  > ° 


This  together  with  the  simple  zeros  of  K(s)  at  v + ~ yields  the  particular  sign 


pattern  iii)  of  (2.10). 


When  t = — (2.9)  be 

2 


comes  a formula  for  f'j-*-  : it  is,  after  multiplication  by 


-1,  precisely  the  formula  given  by  Schoenberg  in  [h)  and  again  in  [9] , as  is  seen  when 

(2.10)  is  evaluated  for  t = More  generally  the  formulas  of  Theorem  2 reduce  to 

formulas  of  C.  de  Boor  and  1.  J.  Schornbcrq  (1)  when 

j even  anti  t = ^ 

or  when 

j odd  and  t - 0 . 

In  the  case  i = 0 and  t = the  formula  has  also  Iren  established  by  C.  A.  Micchelli  15]. 
For  n = 4 and  j - 1 the  formula  is 

(2.9),  E’(t+  t > f»  (t)  - P^(t+|)f’  (t)  = l A^ffv)  + / K(s)f <4)  (c)ds 


-0- 


r 


with  A^  and  K.(s)  given  just  as  in  (2.10)  except  that  ii)  is  replaced 


(2-10) , 


a.  (t)  = 3U 

1 


X - 1 1 + X 

12  l 1 

— t <4t  +1)  + (-1)  5 — — 4t 

\ 1 


i = 1,2 


Ap  = AQ(t)  = 6Vi(4t2+l),  0 < t < | • 

Clearly  A (t)  >0  as  are  a.(t),  0 < t < Thus  the  A of  (2.9)  have  the  desired 

0 i — — 2 v • 

v 

sign  pattern  (-1)  A >0. 

And  as  a last  example  for  n = 4 and  j = 2 we  lave 

00  oo 

(2.9)  K(t+i)f " ' (t)  - E'“  (t+i)f"  (t)  = T A f(v)  + / K(s)f <4)  (s)ds 

3 2 3 2 u v 

— 00  —00 

where  now  ii)  of  (2.10)  is  replaced  by 


1 - X, 


*l(t)  =a2(t) 


-24y 


(2.10), 


1 


AQ  = 48p 


Having  thus  concluded  our  examples  of  some  of  the  formulas  contained  in  Theorem  2, 
we  now  use  the  general  formula  to  prove  Theorem  1.  The  argument  is  along  lines  very 
like  those  used  above  to  derive  the  Proposition  concerning  A from  the  Polya-Szego 
formula  (2.1). 

Proof  of  Theorem  1 . Each  A is  a convex  set  and  so  can  he  completely  described  in 
terms  of  its  lines  of  support.  Just  as  in  the  proof  of  the  Proposition,  determining  the 
position  of  the  supporting  lines  in  a given  direction  with  normal  (a,B)  amounts  to 
maximizing 


(2.11) 


fc’<j)  (t+i)f (j+1)  (t)  - (t+i  )f <j)  (t) 

n-1  2 n-1  2 


over  all  f 6 F . For  f c F we  evaluate  (2.11)  via  (2.5)  as 
n n 


(2.12) 


l A f(v)  + / K(s)f  (n)  (s)us  < l |A  | 4 y / | K (s ) |ds 


where  the  inequality  follows  from  conditions  (1.1)  defining  the  class  F . Clearly 
equality  occurs  in  (2.12)  if  and  only  if  f satisfies  both 


-9- 


f (v)  = signum  = (-1) 


f (s)  = Y signum  K(s)  a.e.  . 
nn 

But  by  comparing  (1.3)  with  (2.6),  we  sec  that  these  last  two  conditions  are 

satisfied  by  the  Euler  spline  E , and  in  fact  these  conditions  characterize  E . So 
J ' n n 

(2.11)  is  maximized  when  f (s)  - E (s)  . This  implies  that  the  pair  (F^'(t),  ''(t)) 

n n n 

is  on  the  boundary  of  A ^ , and  as  we  vary  t we  generate  the  entire  boundary.  We 


mention  in  passing  that  the  ratio 

(t+ 

a n-1 2 

6 = " 'B(i)  (t+f) 

n-1  2 

takes  on  every  extended  real  value  and  so  maximizing  (2.11)  gives  supporting  lines  in 
every  possible  direction.  /// 

From  the  uniqueness  comments  made  above  in  the  proof  of  Theorem  1,  we  obtain  the 
following  remarkable  property  of  the  Euler  spline: 

Corollary.  Assume  n >_  3 and  f c F . The  equations 

f (j) (t)  = F (j>  (t) 
n 

f(j+1)<t)  = A(j4J,(  t) 
n 

can  occur  simul  tan  ■>  .c-l  y .jt  sninr.  poi  nt  t onl  y i f 


f(s)  = F (s) 
n 


for  all  ml  s.  lf  i 0 wo  exclude  from  our  asscr t i on  any  integra  1 va  1 u e of  t. 

In  other  words,  the  pair  (f(i)(t),  f<j0)(t))  is  always  in  the  interior  of  A 
unless  f is  the.1  Euler  spl  i no  t in  which  case  the  pair  is  a 1 ways  on  the  boundary  of 


-10 


13.  A construction  of  the  formulas  of  Theorem  2. 


We  will  carry  out  the  construction  for  n even;  for  the  case  of  n odd,  small 
variations  are  necessary.  Our  main  task  is  to  construct  K(x)  with  the  properties 
given  by  (2. A).  The  formula  (2.5)  then  emerges  easily  by  integration  by  parts. 

There  are  two  main  tools  involved  in  the  construction,  tools  from  the  theory  of 
cardinal  spline  functions.  Our  references  for  this  material  are  |9)  of  Schoenberg  and 
(1!  of  de  Boor -Schoenberg  to  which  we  refer  the  reader  for  details;  also  to  (51  and  [f>] 
where  C.  A.  Micchelli  has  developed  some  of  these  methods  to  provide  an  "optimal 
estimator"  for  f'(t). 

The  eigensplinos.  These  are  cardinal  splines  S satisfying  the  functional  equation 

S (x+1 ) = AS  (x)  . 


The  number  A is  called  the  eigenvalue.  Wo  need  two  classes  of  such  eigenspl ines , 
those  vanishing  at  the  integers  v and  also  those  which  vanish  at  the  points  v + 
in  both  cases  the  knots  are  to  be  at  the  integers . 

When  n = 2m  the  degree  of  K is  2m  - 1 and  according  to  (2.0)  K must  have 
sign  changes  at  v t ^ for  every  integer  v.  We  find  in  [9]  2m  - 1 eigenvalues 

P < ...  < P , < P = -1  < p < ...  < p_  , < 0 
) m-1  m mil  2m- 1 


and  corresponding  eigenspl  ines 


(3.1) 


S . ( x ) , i = 1 , . . . , 2m  - 1 

l 


of  degree  2m  - 1 


(3.2) 


sat i s f vi ng 


S.(x+1)  = p.S.(x)  for  all  x . 


We  note  that  for  i m + 1 , . . . , 2m  - 1 

1 im  S . ( x ) = 0 ; 

l 

x 

and  for  this  reason  these  S.(x),  (i  = nH  1 , . . . , 2m-l ) are  sometimes  called  the  "decreasing" 

l 

eigensplinos.  Among  the  eigen:  pi ines  (3.1),  S (x)  is  the  only  one  which  is  bounded. 


When  n = 2m  -t  1,  we  again  find  19)  2m  - 1 eigenvalues 

X < ...  < X < X = -1  < X < ...  < X„  < 0 
1 m-1  m m+1  2m- 1 


and  corresponding  cigensplincs 
(3.1)  , 

of  degree  2m  satisfying 


S . (x)  , i = 1 , ....  2m  - 1 
1 


(3.2) , 


S. (0)  =0 

1 

S.(x+1)  = X.S.(x)  for  all  x . 

1 li 


The  Budan-Fourier  Theorem  for  spl ines ■ For  a given  spline  function  f,  we  let  Z (a,b) 

denote  the  number  of  zeros  of  f,  counting  multiplicity,  on  the  open  interval  (a,b). 

If  f is  of  degree  n,  f^  is  piecewise  constant  and  Z . . (a,b)  is  defined  as  the 

f(n) 

number  of  strong  sign  changes  on  (a,b);  thus  an  interval  where  f vanishes 


identically  is  ignored.  In  addition 


.(n) 


(n) 


S (f (a) , ...,f  (a) ) 

denotes  the  number  of  sign  changes  in  the  sequence  f (a ) , . . . , f ' ' ' (a ) where  zeros  are 
ignored.  Similarly, 

S+(f  (a) f <n)  (a)  ) 

counts  tho  sign  changes  with  zeros  taken  positive  or  negative  so  as  to  maximize  the  count. 
With  these  notations  we  state  the  useful 

Theorem.  Assume  that  t he  spl  i no  f is.  erf  precise  degree  n and  has  a finite  number 
of  simple  knots  in  (a,b) . Then 

(3.3)  Z,(a,b)  < Z , . (a,b)  + S_  (f  (a)  , . . . , f (n)  (a+) ) - S+  (f  (1.)  , . . . , f (n>  (b- ) ) . 

f - f(n) 


There  are  many  references  to  this  result;  perhaps  the  most  accessible  for  the  present 
purposes  is  (1)  or  (6). 

The  use  of  the  Budan-Fourier  Theorem  in  the  presence  of  cigensplincs  is  very  much 
facilitated  by  the  following  proposition  which  plays  a very  important  role  in  our 
construct  ion. 
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Projaog i i ion.  1.  'Die  eiqonspl ines  s (x)  (3,1)  gatisfy  fog  every  integer 

1 


s"(s. (v*  ~ 


( v-t  2 ) ) = 1 


(3.4) 


S^S.  (v+  * ) s|?m-l)(v+  | )}  = . 


2.  The  eigcnsp]  ines  S ^ (x ) o^f  (3.]),  satisfy  for  every  integer  v 


(3.4) , 


s'(s.  (v),...,S.(2m)  <v+))  i 

1 1 

S+(S.  (v>  , . . .,F(2n°  (v-))  — i . 
1 1 


The  proposition  appears  in  [1]  and  also  (6);  it  is  proved  on  the*  basis  of  the 
Gantmachai -Kre i n Theorem  on  oscillation  matrices. 

Determining  the  kernel  K (x ) . Set  n 2m,  in  _>  2 . Fix  j and  t as  in  Theorem  2.  For 
simplicity  we  assume  0 < t < ~;  clearly  this  represents  no  essential  restriction.  Also 
if  j = 0 we  exclude  t = 0 as  indicated  in  the  remarks  following  the  statement  of 
Theorem  2.  I’ut 

, is  ^ V*  r.  / \ 2m- 1 , , . 2m-l-j  , 2m-2-j 

K.  (x)  - ) a.5.(x)  ♦ ax  4 b(x-t)  J + c(x-t)  J,  x < 1 

1 .c,  1 1 4 4 + — 

i = l 


(3.5)  K (x) 


2m- 1 


K.  (x)  y a.s.(x),  x > t 

V 2 . L , l l ~ 

>-  i s=m  + 1 

for  an  appropriate,  choice  of  the  2m  4 l parameters 


{ a. 


in- 1 m4 1 


2m- 1 


, a , b , r ) 


to  be  determined  presently.  Note  that  in  (3.5)  each  S . (x) , i = - 1,  is  extended 

from  the  interval  (-1,0)  to  (-1,1)  without  a knot  at  0;  instead  the  term  ax^m4* 
provides  the  knot  at  0 for  K(x).  To  check  that  K(x)  is  well  defined  by  (3.5)  as  a 
single  valued  function,  both  definitions  of  K(x)  given  by  K^(x)  and  K^fx)  must  agree 
on  the  overlap  (t,l).  So  when  restricted  to  the  interval  (t.,1),  (x)  and  K^(x) 

must  bo  identically  the  same  polynomial . Equivalently 

Kl<{>  = K2C)  1 = 0 2m  - 1 

for  any  fixed  t witli  t < t < 1.  After  a little  rearrangement,  these  conditions  yield 


a linear  system  of  2ra  equations  in  the  2m  unknowns 


{; 


,b,c} 


with  a right  hand  side  given  by  the  term  ax2m  1 evaluated  at  x = t. 

To  determine  a solution  of  the  above  linear  system,  consider  first  the 
homogeneous  system  obtained  by  setting  a = 0.  Suppose  there  were  a nontrivial  solution. 
Then  the  result  is  a K(x)  defined  by  (3.5)  on  the  entire  axis  but  with  no  knot  at  0 
since  a = 0.  Now  it  is  easily  seen  that  the  sum 


2m- 1 


(3.6) 


K(x)=  ) a.S.(x),  x > t 

i=iiHl  1 1 


is  nontrivial.  In  addition,  the  functional  equations  (3.2)  and  the  ordering  of  the  eigen- 
values ii.  together  imply  that  for  large  values  of  the  argument  x the  sum  (3.6)  is 

dominated  by  S (x)  . Thus  from  (3.4)  witli  v large  we  have 
m 1 1 

(3.7)  S+(K(v+|)  ,K'  (v+ j)  K2m_1(v+^))  < m * 1 . 

Similarly 

(3-8)  s"(K(-v+  * ),K'  (-V+-*  ) K2m_1(-v+i))  < m - 2 . 

Now  using  these  two  estimates,  we  apply  on  each  of  the  intervals  ( - v+ - , t ) and 
(t,v+  - ) the  Budan-Four ier  Theorem  (3.3)  to  K(x)  given  by  (3.6)  with  a = 0.  When 
the  resulting  two  inequalities  are  added  together,  we  obtain 

2v  - 1 < 2 V - 1 + s‘(K(tt>  ,K'  (t-O K<2m-1  ) ((_f)  , _ s+  (K  (t-)  ,K'  (t-)  , . . . (t-)  ) 

(3.9)  + (m-2)  - (m+1) 

< 2v  - 1 + 2 t (m-2)  - (ml)  = 2v  - 2 . 

The  second  inequality  of  (.1.9)  follows  because  the  two  sequences 

K ( t+) , K ' <t+> K(Zm"1) (t+) 

K(t-)  , K ' (t-) ,K<2,n'1’  (t-) 


can  differ  (by  (3.5))  only  in  two  consecutive  entries:  hence  the  corresponding  difference 
in  (3.‘>)  is  at  most  2.  Now  (3.9)  is  a contradiction,  implying  that  the  homogeneous 
system  has  only  the  trivial  solution. 
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Now  set  a = 1 and  so  obtain  a unique  K(x)  defined  by  (3.T>).  To  this  function 
K(x)  wo  again  apply  the  above  arguments  leading  to  (3.9)  but  now  with  the  one  change 
that  0 is  a knot.  The  result,  valid  for  all  large  integers  v,  is 

(3.10)  2v  - 1 < 2v  + 2 t (m-2)  - (m+1 ) 2v  - 1 . 

So  we  must  have  equality  in  (3.10)  which  forces  equality  in  (3.7)  and  (3.8).  From  these 

equalities  we  can  easily  derive  all  t.ho  properties  asserted  for  K(x)  in  Theorem  2. 
Properties  of  K(x).  ii)  of  (2.6)  is  clear  from  (3.5),  as  is  the  exponential  decay  of 
K(x).  From  (3.2)  and  (3.'>),  K(vt  ’ ) = 0 for  all  integers  v;  that  these  zeros  are 
simple,  and  that  K(x)  has  no  other  zeros,  follows  from  (3.10).  Thus  K(x)  changes 
sign  at  each  point  v 4 * . 

Again  from  (3.10),  K(?m~1)(x)  must  change  sign  across  every  integer  and  when 

j > 1 these  are  clearly  t lie  only  sign  changes  of  K ^ * (x)  . For  j = 0 there  is  a 

possible  sign  ch mgr  at  t,  but  we  will  eliminate  this  possibility  shortly.  Foimula 
(2.r>)  emerges  by  integrating  by  parts 

J K (x)  f ^ (xhlx  . 


T)i us  the  A of  (,’.r>)  are  given  by 
v 


(3.11) 

So 

(3.12) 

and  we  nomtalizr  by 


, <?m-l ) , , „ ( 2m- 1 ) . , . 

A - -(K  (v+)  - K (v-) ) 


A strictly  alternates  in  sign 
v 


sign  A^  > 0 . 


This  normalization  imp  I !•  •' 


K ( 2tn- 1 ) (x)  > 0 for  .j  < x < o 


(x)  > (l  for  -2v  - 1 < x < -2v  ; 


in  part icular 


K^Wi,  > 0 . 


(3.13) 


Now  equality  in  (1.8)  combined  with  (1.12)  yields 


(1.14) 


sign  K ' (-2v-  ~ ) = (-l)m  2 = (-l)1"  . 


(-l)mK(x)  > 0 for  -2v 


1 

— < 
2 


x < -2 v + 


1 

2 


and  due  to  the  simplic'  y of  the  zeros  of  K(x)  we  find 
(1.15) 

valid  for  all  v.  Th. s - 'tablishes  iii)  of  (2.6). 


-l)m4VK(x)  > 0 for  v-~<x<v  + ^ 


Concerning  the  case  j = 0,  we  sen  from  (1.5)  that  K has  a double  knot  at  t; 
this  allows  a possible  change  of  sign  in  K^2m  at  t,  and  we  must  eliminate  this 
possibility  in  order  to  preserve  (3.12).  Given  the  sign  changes  of  ^ at  every 

integer,  a sign  change  at  t would  entail 

, ( 2m- 1 ) 


K 


(2v-i)  < 0 


for  large  positive  v.  Following  the  same  line  of  reasoning  which  resulted  in  (3.15), 
we  would  arrive  at 

(-l)m,VK(x)  < 0 for  v— ^<x<vt~. 

This  contradicts  (1.15);  hence  there  is  no  sign  change  of  K*‘m  at  t. 

Thus  we  have  established  formula  (2.5)  with  a right  hand  side  described  by  (2.6). 
From  (1.5)  and  the  integration  by  parts,  it  is  clear  that  the  left  hand  side  of  our 
formula  is  of  the  form 

of  (:i)  (t)  + Rf  ^ + 1 * (t)  . 

We  have  yet  to  determine  a and  6,  or  more  precisely  the  ratio  R/a,  as  our 

formula  is  determined  only  up  to  a multiplicative  constant. 

Recall  the  sets  A.  of  51.  For  every  s 
3 

(f;(j)  (s),£h+1)  («) > * a.  ■. 
n n j 

and  in  fact  on  the  basis  of  all  the  properties  (2.6)  of  formula  (2.5)  and  the  correspond- 
ing properties  (1.3)  of  E (s) , we  can  already  conclude  as  in  12  that 


-16- 


Wi  th 


(£'(j)  (n)  (r.)  ) f ZA 

n i) 


= E 


,(j) 


(s)  , 


y =£(j+1),s) 
n 


we  f i nd 


dy 

^ = dT 

dx 


j / :ix 
/ ds 


ff(i+2)(s) 
n_ 

F(^41) (s) 
n 


n-I 2 _ 

P<;1!  (S4  i) 
n- J 2 


And  from  Figure  1 it  is  clear  that 


dy 

dx 


:=t 


So  we  have  0 = £(j)  ( t+ ~ ) and  a = -£'(:’+1)  (t+  ) . 

n-l  2 n-l  2 

The  odd  case  n = 2m  + 1 is  settled  in  exactly  the  same  way  with  the  even  degree 
eigensplines  (x)  given  by  (3.1),  and  (3.2),  replacing  the  (x)  . One  then  argues 
on  the  integer  points  x = v,  as  indicated  by  (3.4),« 
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